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Abstract 

The various inadequacies of Callaway's phenomenological 
model of lattice thermal conductivity has been critically 
analyzed and the model is repaired in a modified form in 
which the systematic replacement of life time by line widths 
amicably resolves the various issues. The involvement of 
various scattering events in the heat transport, e.g., 
boundary scattering, impurity scattering, anharmonic 
phonon scattering, resonance scattering and interference 
scattering has been addressed in the new framework with 
the help of quantum dynamical many body theory. The 
technological importance of Ge is well known and hence it 
becomes significantto investigate the thermal behavior of it 
in details as its electrical properties are temperature 
dependent. Further, the CdTe also shows its vital 
importance in the fabrication of infrared optical windows 
and photo voltaic solar cells. The phonon heat conductivity 
of Ge and CdTe in the temperature range 3.3-298 K and 
1.780-239.260 K based on the modified Callaway model have 
been analyzed and excellent agreements between theory and 
experiments are reported. The present formulation is found 
to be well justified and can be successfully applied for the 
calculations of thermal conductivity of several other 
crystalline solids. 
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Introduction 

Aconsiderable interest throughout the globe for more 
than a decade has been witnessed to develop the 
semiconductor based optoelectronic and nano 
technological devices. Of the widely used 
semiconductors Ge and Si, their unique physical 
properties like indirect band gap and large effective 
masses, the effort to develop mid infrared and fiber 
optic communication devices was limited (Wang K L, 
Cha D., 2007}, but some efforts have been pioneered to 
design low cost Ge-based photodetectors for near far 
infrared communication bands and high resolution 
spectral imaging (Wang J, Lee S., 2011).The modern 



low dimensional microelectronic devices on the same 
substrate such as quantum wells in the Si/Ge systems 
and Si/Ge and Si/Ge/C quantum dots have been 
tailored recently (Egorov V A, Cirlin G.. 2004). On the 
other hand, direct optical band gap (~1.5eV), III- VI 
compound semiconductor CdTe has been recognized 
as a very promising photovoltaic material to fabricate 
infrared optical windows (waveguides) (Rams J., 
Sochinskii N V., 2000) and solar cells (X Wu., 2004) 
and several applications has been developed in the 
field of nanotechnology as artificially tailored 
materials (Novoselov K S.and Castro A H., 2012). In 
semiconductors, the lattice thermal conductivity has 
been emerged as a fundamental transport property to 
facilitate the understanding of heat dissipation in the 
field of nanoelectronic and optoelectronic device 
fabrication with increased efficiency. 

In the crystalline solids where most of the heat is 
transported by phonons, the theory of lattice thermal 
conductivity based on the Boltzmann transport 
equation approach via relaxation time approximation 
has sufferedfrom the commonshortcomings of kinetic 
theory (Peierls R., 1997). A more rigorous basis was 
acquired by the phonon transport theorywith the use 
of energy flux correlation techniques (Kubo R., 1957), 
however,these theories could not probe into successful 
numerical estimates. Of the various theories Callaway 
presented the most successful numerically amenable 
model to analyse thermal conductivity of non metallic 
solids. Following Callaway's theory, the expression for 
lattice thermal conductivity is given by (Callaway J., 
1959). 

K = {k B Tll7i\ p )(pK) 2 f T(a))a) 4 e ma> (e mo) -l)~ 2 do) (1) 

Where v p ,co D ,T{co) stand for phonon velocity, Debye 
frequency, total relaxation time respectively and 
P = (k B T)~ l .In this model, Callaway made some 
assumptions, namely; (a) the phonon dispersion and 
anisotropy effects were ignored, thus the longitudinal 
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and transverse phonon modes behave identically, (b) 
the relaxation time has been considered as a function 
of frequency and temperature except boundary 
scattering relaxation time and (c) the relaxation times 
for different scattering mechanisms have been 

inversely added as T 1 (&>) = X 1 , where T; 

i 

represents the relaxation time of the process Y. This 
greatly oversimplified model is, however, not 
complete and rigorous and leaves immense 
possibilities for modifications via above cited 
assumptions. In this context, Holland(Holland M G., 
1963) made pioneering efforts in resolving Callaway 
expression of thermal conductivity into longitudinal 
and transverse components which werefurther refined 
by SDV model (Sharma P C, Dubey K S, and Verma G 
S., 1971), Tiwari model (Tiwari M D, and Agrawal B K., 
1971) and others (Gairola R P., 1983, Bahuguna B P, 
Painuli C P and Indu B D.,1991). SDV and Tiwari- 
Agrawal models tried to repair the Holland model by 
including some generalized dispersion relations and 
found excellent agreements with experimental 
observations. Since the relaxation time t) 

appearingin the heart of the problem has been worked 
out in many ways by several physicists (Pomeranchuk 
I., 1941, Klemens P G., 1965, Klemens P G., 1958, Erdos 
P, and Haley S B 1969) for a very long time, but in 
most of the treatments the frequent violation of 
Matthiessen's rule has beenobserved. In the modern 
approach (Gairola R P., 1983, Bahuguna B P, Painuli C 
P and Indu B D., 1991, Ansari M Ataullah, Ashokan V, 
Indu B D., 2012) the problem of evaluation of 
relaxation times has been taken up with the new 
perspectives in which the problem of exclusion of 
dispersion and violation of Matthiessen's rule has been 
amicably resolved. 

In the present paper, these inadequacies have been 
overcomewith the help of quantum dynamical phonon 
Green's function theory (Gairola R P., 1983, Bahuguna 
B P, Painuli C P and Indu B D., 1991) in which the 
relaxation time hasbeen evaluated from phonon life 

time. On the basis of the new model, the thermal 
conductivities of Ge and CdTe has been analyzed and 
excellent agreements between theory and 
experimental data have been obtained in the following 
sections. 

Theory 

Owing to Matthiessen's rule, the relaxation times for 
various scattering events, namely; boundary scattering 



(t b ) , impurity scattering (t d ) , cubic and quartic 
phonon scattering (^L.^L) dislocation scattering 

( T dis~) ' resonance scattering (r^ 1 ) , etc. can only be 
considered inversely additive if one considers them as 
independent (non-interacting) processes. However, it 
is not possible to consider their non interacting 
existence, because, there is sufficient possibility of 
interaction of these different collision processes. 
Contrary to this reality, we have come across to the 
most of the work which frequently violates the 
Matthiessen's rule. This problem can be resolved by 
making use of the relation between relaxation time T 
and phonon frequency time width T k (co) as (Gairola R 

P., 1983, Bahuguna B P, Painuli C P and Indu B D., 
1991, Ansari M Ataullah, Ashokan V, Indu B D., 2012, 
Indu B D., 1992). 

T~\(0) « r k (0)) + T c l B +Tft (2) 

Where the phonon line width T k (co) can be written as 
a sum of individual contributions to phonon line 
width coming from defects T®(a>) , Cubic 

anharmonicities T k A (<o) , quartic anharmonicities 
T 4 k A (a>) and impurity-anharmonocity interactions 
T^ D (a>) in the form 

r t (®) = rf (®) + rf (®) + r; A m + r/ D («) + T £ + t r ] (3) 

This treatment suggests the two fold modifications in 
Callaway model, namely; (i) the problem of violation 
of Matthiessen's rule is resolved through phonon line 
width concept, and (ii) T k (w) itself describes the 

details of frequency spectrum hence the inadequacy of 
exclusion of phonon dispersion has been overcome. 

The boundary scattering relaxation time t b = v p I L 

( L = Casimir's length, which has been used as a 
parameter without any physical justification) can be 
modified by considering the involvement of crystal 
microboundaries and micro scale fluctuations in the 

form Tg = v I L(B) where L(B) the modified 

Casimir's length includes the effects of 
microboundaries(Indu B D., 1980, Ansari M Ataullah, 
Singh Nitin P, Indu B D., 2007). 

The phonon line width (co) has been evaluated with 

the help of double time thermodynamic Green's 
function theory via an almost complete Hamiltonian 
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which comprises the effects of defects, anharmonicities 
(Gairola R P., 1983, Bahuguna B P, Painuli C P and 
Indu B D., 1991, Ansari M Ataullah, Ashokan V, Indu 
B D., 2012, Indu B D., 1992) and the various 
components can given by 

T k (co) = 8^(o)X R(~k,k l )R* (-k,k l )co ki S(a) 2 - a>\ ) 



~ AjO 4 + A 2 co 2 (4) 



Here A l and describe the mass change and force 

constant parameters, respectively which provide 
additional information as compared to Klemen's 
results containing the effects of mass change only 
(Klemens P G.,1958). The importance of this new term 
cannot be easily ignored because a very small change 
in central or noncentral force constants drastically 
changes the specific heat of a crystal [note that the 
specific heat term is an inevitable component of lattice 
thermal conductivity expression; 

- { v 2 cos 2 (<?) a (k)C ph (k)d 3 k . 



K 



{in) 



Since the number of excited phonons continuously 
increases with gradually rising temperature towards 
LTC maximum, the probability of phonon-phonon 
scattering becomes comparable to defect scattering.If 
adequate information about the specimen as required 
in the first part of this equation is available, one can 
make parameter free calculations specially in case of 
tailored materials. Also 



R(k v k 2 ) = (w k Jw ki )c(k l ,k 2 ) + D(k l ,k 2 ) 
+ 41 C , k 2 )d (-k 2 , k 2 ) co k 



(5) 



The study of the effects of phonon-phonon and quartic 
phonon scattering has been a matter of long debate 
and the frequency and temperature dependent nature 
of these processes have been varied at the will of the 
workers and remained unjustified(Callaway J., 1959, 
Holland M G., 1963, Pomeranchuk I., 1941, Klemns P 
G, 1965, Klemns P G, 1958, Erdos P, and Haley S B., 
1969). In the present work, this has been sorted 
out(Gairola R P., 1983, Bahuguna B P, Painuli C P and 
Indu B D., 1991, Ansari M Ataullah, Ashokan V, Indu 
B D., 2012, Indu B D., 1992) to give 



t\ a (co) = XVco k 6(co L -co k )l \6xN ptia Q v" * BaT (6) 



T 4 k A (co) = {Til A%M){aJiri(p IV V / 4tt 2 j6<p !I v 3 f 
x(a> 2 k + a\ - a> 2 k+kx ) « B u co 2 T 2 



(7) 



where A is the dimensionless quantity and 

= co L sin(;rfca ) which shows the exact frequency 

and temperature dependence. The various symbols 
appearing in the above equations are obtainable in the 
form (Indu B D., 1992) and 

= % ± \ (8) 

S ±j3= 1±n k l n k 2 +n k 2 n k 3 ±n k 3 n k l W 
a>±a = ®h ± ®h ' a ±P = ®h ± ®h ± ®h (10) 



7,-- 



®k®k 2 ®k, 



cb k ^Pha k 
m =— coth- 



2 (11) 



All of the above results are based on non-perturbative 
approach. 

When the phonons of anharmonic fields enter the 
localized fields and interact with the phonons of 
impurity field, the impurity anharmonicity interaction 
modes occur which can be recognized by the line 
widths 



rf (a>) = {UVh_ 2 ItiP h a )[(M 2 c(l-c)/4NjuV) 2 ]co 4 k 

Yf (a) = (3h/j_ 2 1 MN)[(Mla h<p' V V 1 16a- 2 Pep 1 ' 'v n) 2 
xc(l-c)]a> k * D'co 4 T 2 



(12) 



(13) 



h = t cb + A i a>A + A 2 ffl2 +( B + B i e~ 0laT )m 2 T + Dm 4 T + t r 1 (14) 

It is wise to note here that the processes related to 

quartic anharmonic interactions, e.g., F k A (<y), T k D {co) 

contribute tohigh temperatures and will not be 
included in the present study. 

The presence of bump or dip in the thermal 
conductivity curve can be explained by a typical 
resonance phenomenoncharacterized by the relaxation 
time (Pohl R O., 1962, Walker C T, and Pohl R O., 1963) 



Rco 2 T n 



(15) 



(cvl-co 2 ) 2 +{CLlnfco 2 cvl 
where R is a proportionality constant containing the 
concentration of impurities causing the resonant 

scattreing, &> is the resonance frequency, and Q 
describes damping of the resonance. 

Analysis of Thermal Conductivity 

For the purpose of analysis, we have taken the 
experimental values of the work of Hamilton-Parrott 
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and Holland (Hamilton RAH, and Parrott J E., 1969, 
Holland M G., 1964). For the sample of Ge, the 
observations are available in the temperature range 
(3.3-298K) and those for CdTe are reported in the 
range (1.780-239.260K). Due to the lack of sufficient 
information about the samples, we have taken some of 
the quantities in the analysis as reasonably 
justified(Gairola R P., 1983, Bahuguna B P, Painuli C P 
and Indu B D., 1991) parameters which are furnished 
as Table [1]. It is observed that at low temperatures, 
below the maximum thermal conductivity due to the 
excitation of low frequency phonons, the major 
contribution to thermal conductivity comes from the 
combined boundary scattering. As the temperature 
starts readily to rise, phonons of relatively higher 
frequencies get excited and start interacting with point 
impurities present in the specimen giving rise to 
defecting scattering phenomenon. Owing to gradual 
rise in temperature nearmaximum of the thermal 
conductivity curve the higher frequency phonons 
become available, which, in turn, invokes the phonon- 
phonon scattering processes. At the same time the 
phonons present in the anharmonic fields start 
interacting with localized phonon fields of point 
impurities giving rise to impurity-anharmonicity 
interaction modes. Thus the phonon-phonon 
scattering processes and interference scattering 
dominate in the close vicinity of the thermal 
conductivitymaximum and above. 

Electrically inactive monoatomic impurities also give 
rise to the resonance scattering in the close vicinity of 

(Oq which are also found to be very important in the 

present samples and capable to provide the best fit at 
the dip in the conductivity curve. The present 
calculations depicted in the figures [1] and [2] show an 
excellent agreement with experiments. 

TABLE 1 CONSTANTS AND PARAMETERS 



Sample 



L(B) [cm] 



A xlCT 44 [S 3 ] 



A 2 xl0 23 [S] 



BxlQ 23 [SK '] 



B, xlO 23 [SK l ] 



DxlO U [S 3 K '] 



flxlO [c/s] 



a> Q xlO [c/s] 



Q [c/s] 



a[cm] 



Ge 



0.146429 



0.022 



20.734 



499.09911 



53.85572 



1.062 



1.421 



6.05551 



1.79 



CdTe 



0.0861 



0.110 



598.09911 



22.998 



5.32350 



0.41 



v[cms ] 


3.11xl0 5 


3.11xl0 5 


D [K] 


118 


252 




Temperature ( K) 

FIG. 1 ANALYSIS OF THERMAL CONDUCTIVITY OF Ge 
SAMPLE 




50 100 200 



Temperature ( K) 

FIG. 2 ANALYSIS OF THERMAL CONDUCTIVITY OF CdTe 
SAMPLE 

Conclusions 

It can be statedfrom the present study that the various 
inadequacies left in the Callaway model havebeen 
almost removed by the quantum mechanical treatment 
of the phonon line widths. This modified model 
successfully explains the experimental observations of 
thermal conductivity and can be applied for the 
calculations of thermal conductivity of several other 
optoelectronic materials including artificially tailored 
low dimensional solids. 
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